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Dark energy and dark matter are only indirectly measured via their gravitational effects. It is 
possible that there is an exchange of energy within the dark sector, and this offers an interesting 
alternative approach to the coincidence problem. We consider two broad classes of interacting 
models where the energy exchange is a linear combination of the dark sector densities. The first 
class has been previously investigated, but we define new variables and find a new exact solution, 
which allows for a more direct, transparent and comprehensive analysis. The second class has not 
been investigated in general form before. We give general conditions on the parameters in both 
classes to avoid unphysical behavior (such as negative energy densities). 

PACS numbers: 98.80.-k 95.35,+d 95.36,+x 98.80.Jk 
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I. INTRODUCTION 

Cosmological observations point to the existence of 
non-baryonic cold matter and of a late-time acceleration 
of the universe (see, e.g. d 0, [|). If gravity is modelled 
on cosmological scales by general relativity, and if we as- 
sume that the universe is homogeneous and isotropic on 
these scales, then the late-time acceleration is sourced by 
a dark energy component, and the universe is dominated 
by the "dark sector" . The quest to uncover the true na- 
ture of dark matter (DM) and dark energy (DE) is one 
of the most pressing topics of modern cosmology. 

One of the fundamental puzzles within this quest is 
the "coincidence problem": how is it that we seem to 
live in a time when the densities of DM and DE are of 
the same order of magnitude, given that they evolve very 
differently with redshift? An interesting proposal is that 
interaction between the dark fields could perhaps allevi- 
ate the coincidence problem. Various interaction models 
have been put forward and studied (see, e.g. Hi p a, [ a EL la 

im m m m m he he m m m he ^i^iiiiii). 

For a flat FRW universe, the background dynamics 
after recombination are governed by the equations of en- 
ergy balance and the Raychaudhuri field equation: 

pb = -3Hp b , (la) 

Pc = -3Hp c + Q, (lb) 

p x = -3(1 + w x )Hp x - Q , (lc) 

H = -AwG [p b + p c + (1 + w x )p x ] , (Id) 

where H = a/a is the Hubble parameter, p c is the cold 
DM density, pb is the baryonic density, p x is the den- 
sity of DE and w x < is its constant equation of state. 
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The baryons only interact gravitationally with the dark 
sector, and Q is the rate of energy transfer in the dark 
sector. The Friedmann constraint equation is 



2 8irG . . 

H = — r— (Pb + Pc + Px) 



(2) 



Note that the field equations (|ld[) and ^ are indepen- 
dent of Q, because of total energy conservation. A posi- 
tive Q > represents transfer of energy from DE to DM; 
a negative Q < represents transfer of energy from DM 
to DE. 

In this paper, we consider interactions that are linear 
combinations of the dark sector densities: 



Q = A c p c + A x p x 



(3) 



Here the rate factors Ai are either proportional to H or 
constants, leading to two classes of interaction model: 



Model I Ai = 3atiH , 
Model II At = 3Ti, 



(4) 
(5) 



where ai are dimensionless constants and Tj are con- 
stant transfer rates. Observations impose the general 
constraint that the interaction should be sub-dominant 
today, so that 



H<i, |r/|< So- 



ft) 



Model I has been recently analysed by Hi and earlier 
work considered the special cases a c = a x [|| and a x = 
We use an alternative approach, defining new 

variables to simplify the parameter space, and finding the 
general exact solution. We are able to recover previous 
results more directly and simply, and to provide some 
new insights into the model. 

The H in the Q-term for Model I is motivated purely 
by mathematical simplicity. By contrast, the energy ex- 
change in Model II is motivated by similar models that 
have been used in reheating (25[, curvaton decay [26[ 
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and decay of DM to radiation [27fl . As far as we know, 
Model II for the dark sector interaction has has not 
been treated in the general case before. The special case 
Tj; = has been analyzed by Q (and by 0] in the case 
where DE is modelled by exponential quintessence). 

A summary of the paper is as follows. Some properties 
of the general case (Q not specified) are presented in 
Sec. HU Model I is studied in Sec. [TTTJ In Sec. Ml we 
analyze Model II. Finally, we conclude in Sec. [V] 



II. THE CASE OF GENERAL Q 



We define the dimensionless dynamical variables [2S 



8ttG 



8nG 



8ttG 



Q- (7) 



Before proceeding further, we now show that we can 
"hide" the constant DE equation of state by defining a 
new interaction variable z := z/(-Wi); in this way, all 
our results below will be independent of the value of w x . 
The dark sector balance equations read 



y 



3x (1 — x) — z , z := , (8a) 

w x 

-3xy + ~z , (8b) 



where a prime denotes d/d(—w x N), with N — In a (we 
choose ao = 1). The baryonic density is determined by 
the Friedmann constraint, 



8ttG 
3lP Pb 



1 



x-y. 



(9) 



We cannot solve the system of equations until Q is 
specified, but we can draw some conclusions for a general 
Q or z. The simplest cases are z = z(x,y), when the 
system ([5]) is closed and autonomous. Model I is such a 
case. The next simplest cases are those for which z is not 
determined algebraically by x and y, but does satisfy an 
equation of motion of the form z! — F(x,y,z), so that 
we have a 3-dimensional autonomous system. Model II 
is an example of this case. 

The critical points (rr*, y*) of the dynamical system ([5]) 
must comply with the conditions 

3x* (1 — — y») = , (10a) 
— 3x* (1 — + y*) + 2z* = , (10b) 

where z* is the interaction variable z evaluated at the 
critical points. Equation (|10a[) implies that either x* = 
0, which represents the usual matter dominated point, 
or x* + y„ = 1, which implies no contribution from the 
baryonic component. 

From Eq. (fl0b|) . we see that the option x* = di- 
rectly implies that z* = 0. Thus pure matter domina- 
tion can only exist if the interaction term vanishes at 
the corresponding point. On the other hand, the option 
x* + y* = 1 leads to 



If the dynamical system is autonomous, the above equa- 
tion depends only on x* and gives the position of the 
critical point that is compatible with a nonzero interac- 
tion term (and no baryonic contribution). By contrast, 
if the system is not autonomous, we require extra infor- 
mation before we can completely determine its critical 
points. 

Another quantity of interest is the dark energy-to-dark 
matter (DE-DM) ratio, R = p x /p c - From Eqs. dH) we 
obtain the evolution equation 



R' q {x + v) ~ 

— = 3 z . 

R xy 



Px X 

R= — = -. 

Pc V 



(12) 



which is equivalent to Eq. (3) in [291, |30j . As we show be- 
low, this equation can lead to exact solutions in some 
particular models. A key point relates to the factor 
(x + y) in Eq. I|12p : its substitution by the Friedmann 
constraint © in some cases can hide the existence of 
exact solutions (e.g. 0, H^, [IFj, HH ) . The standard non- 
interacting case, z = 0, shows the expected exponential 
result R = R e~ 3w * N , where R = fl x0 /fl c0 - 3.4 is the 
present value of the DE-DM ratio. 

If the current value of R is close to an asymptotic value, 
Rq ~ Rqo := i?(oo), then the coincidence problem is 
"solved". (Strictly, one has transferred the coincidence 
problem to a problem of explaining the dark sector in- 
teraction.) In this case, R must be a slowly-evolving 
function as N — > oo, i.e. R' ~ at late times, which 
requires 



ixy — z(x + y) ~ . 



(13) 



As we shall discuss below, for some models this condition 
is met only at some points of the phase-space. 



III. MODEL I: Q = 3H (a x p x + a c p c ) 

This model was intensively studied in [i[ (our ai corre- 
spond to their A/). We find the critical points in a more 
direct way, and also find the general exact solution. This 
allows us to recover many of their results more directly 
and simply, and to provide some new insights into the 
model. 



A. Critical points and their properties 



The dimensionless interaction variable z is 

a i 



z — 3 (a x x + a c y) , 5/ 



(14) 



3a;* (1 



(11) 



Thus z = z{x, y) and the phase-space is 2-dimensional 
and autonomous. 

As discussed before, the first type of critical point in 
Eqs. (fTO)) is that for which x* — and z* = 0. Such 
a critical point in the present model also needs = 0, 
and then it corresponds to a baryon dominated point. 
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This point is not realistic and we will exclude it from our 
analysis; for more details, see Q. 

The critical points that are compatible with a non-zero 
interaction term comply with the constraint x* + y* = 1 
and are solutions of Eq. (fTTj) in the form 



X ;fc 1 X ^ — CXg^X^ I Q^c fl ^ * J 



(15) 



The critical points and their linear stability are summa- 
rized in Table I; the results are in agreement with those 
of Q ■ For convenience we have defined the parameters 



Ci 

c 2 



- (1 - a x - a c ) , 



C\ - a x a c 



(16a) 
(16b) 



Note that the existence condition for the critical points 
ensures that C 2 is real. 

For a physically viable model, one of the critical points 
should correspond to a DM dominated universe, x — > 
and y — > 1, at early enough times; moreover, this critical 
point should be an unstable point. The only candidate 
is point A, but xa ^ 0. For DM domination we need 
\xa\ <C 1; to linear order in x*, we obtain from Eq. (fT5|) 
that 



x A 



1 - &x 



VA 



2Ci 
1 - a x 



(17) 



At late times we must get a DE dominated universe, 
which should correspond to the stable point B. For the 
latter, Eq. (fH))) gives the estimate 



x B 



2Ci 



VB 



1 - OLr 



(18) 



For certain parameter values the DE density is negative 
at early times [H (see also Q in the cases a x = and 
oe x = otc). In a first approximation we can use Eqs. (fT7|) 
and (p~8j) to determine the conditions under which both 
DM and DE are non-negative and well behaved at all 
times. 

At early times, we must impose the constraint 



< x A 
which is satisfied if 



< 1 . 



<5 C > , 1 > a x + a c 



(19) 



(20) 



or 



a c <0, l<a x +& c . (21) 
Likewise, at late times we have the constraint 



0<y B 



1 — a c 



< 1, 



(22) 



and then the same conditions (|2"D|) and (|2"Tj) apply, but 
now with the interaction constants interchanged, a c <-> 



a x . It can be verified that the conditions described above 
give the correct description for the different cases de- 
picted in Figs. 7, 8 and 9 in |4|. 

It follows that for a plausible scenario, the interaction 
parameters a c and a x must be both non-positive or both 
non-negative. A remarkable result arises now. If the 
interaction parameters are to have the same sign, the 
only permitted case is that of Eq. (|20|) . as the case in 
Eq. (j2Tj) cannot be consistently satisfied if the interaction 
parameters are both negative. 

Another constraint appears from the condition that 
the critical points are such that < x*,y* < 1. From 
Eqs. (fTTl) . we get a x < 1; likewise, from Eqs. JT5|) we 
get a e < 1. Finally, the critical points should be real, 
Cf > 0, and then the allowed values (a c ,a x ) are those 
located below the rotated parabola 



(a x - a c ) 2 > 2{a x + a c ) - 1 . 



(23) 



Therefore, we arrive at the overall conclusion that 
< a c , a x < 1 if the DM and DE contributions are 
to be positive and well behaved, < x, y < 1, at all 
times. This is in agreement with the results of [J] (see 
the critical points and B 2 in Table I). We used a first 
order calculation based on Eqs. (fTT|) and (| 1 8j) ; however, 
by continuity, we expect the result to be true even in 
cases where DE and DM contribute significantly at early 
and late times. 

A more general statement exists regarding the positiv- 
ity of the dark components at early times. Note that we 
can write Eq. (I15[) for the critical points in the form 



(24) 



3w x x*y* = -3 (a x i:, + a c y*) = -z* 



If both DM and DE are to be positive for all times, then 
z and w x should have opposite signs. This explains why 
the DE component becomes negative at early times in 
certain models: it is because the interaction variable z 
can become negative at DM domination. In our case, 
non-negative interaction constants ensure that both DM 
and DE are positive at all times (see also Figs. 7 and 9 
in Q for some other examples.) 

In Fig. [I] we show the allowed values of the interac- 
tion parameters (represented by the red shaded region) 
according to the discussion above. 



B. DE-DM ratio: exact solution 



With the interaction term (|T4j) . Eq. (|T2|) becomes 

R' = -3 [a x R 2 -(l-a x -a c )R+ a c ] . (25) 

This equation provides the known solutions of simpler 
cases, e.g. for the case a x = 0, we have 

R . = /f^o + _^_\ e -S( Wx+ac )N _ (26) 

V ilco w x +a c I w x + a c 



Point 




y* 


Existence 


Eigenvalues 


otability 


A 


G\ + a c — C2 


l — XA 


cf >0 


-3(Ci + a c -C 2 ); 


Unstable if xa < and C2 + a x > 
Saddle if xa > and C2 + Qi > 

or: xa < and C 2 + < 
Stable if xa > and C2 + a x < 


B 


Ci + a c + C 2 


l—XB 


c\ >0 


-3(Ci +S C +C 2 ); 
-6C2 


Saddle if x B < 
Stable if x B > 



TABLE I: Critical points of Model I and their stability. Ci,C2 are given in Eqs. (|16p : notice that C2 is positive by definition. 
The existence conditions are given here in general, but the discussion in the text suggests that both interaction parameters 
should be positive, see Eq. (I20|l . so that both xa and xb are positive too. Thus, point A is of the saddle type whereas point B 
is stable. 



R = 3.4 
FT =0.1 




a x 

FIG. 1: The red shaded region contains the allowed non- 
negative values of a c and ot x satisfying the reality con- 
straint (|23[) . The cyan (dotted) line represents the equality 
in Eq. (|34|l for FU — 0.1, the blue (dashed) line represents 
the equality in Eq. (136 ^ for Rq = 3.4, and the green (dashed- 
dotted) line is Eq. (|37| with Rac, m = 10. The yellow (solid) 
line corresponds to the particular case a c = a x d. The only 
region that provides reasonable values of the DE-DM ratio at 
both early and late times is the one surrounded by the thick 
black (solid) lines. The value Ri — 0.1 was chosen for pre- 
sentation purposes; the allowed region would become much 
smaller for a more realistic value Ri <0.1. 



which directly recovers Eq. (7) of [9(. Here we present 
the new solution in the general case a x 7^ 0: 

R(N) = ^ - — tanh [3w x C 2 (N - N r )} , (27) 
a x a x 

where C\ and Ci are given by Eqs. P^|) . and the integra- 
tion constant is 



Using Eq. (f27|) we can integrate the DM energy balance 
equation to obtain the DM density as a function of N: 



Pcoe 



(-3- 



-d)N 



cosh3C 2 (iV - Nx) 
cosh3C 2 iVi 



1/3 



(29) 



and then the DE density follows directly as p x (N) — 
R(N) Pc {N). 

The asymptotic values of the DE-DM ratio from 
Eq. ([2T]) are directly related to the values inferred from 
the critical points in Table HI 

i?_oo := R(-oo) = Cl ~ C2 = , (30a) 
a x yA 



R(oo) 



Ci +C 2 _ xb_ 

OL x IJB 



(30b) 



It follows that the smallest value of the DE-DM ratio 
R-00 is determined by the critical point A; likewise, its 
largest value i?oo is determined by the critical point B. 
The correspondence between the asymptotic values of R 
and the critical points A and B is not surprising after all, 
because the roots of R' = in Eq. ([25]) are actually the 
ratios inferred from the critical points. 

There is an interesting point concerning the initial val- 
ues of the DE-DM ratio. For given values of the free 
parameters of the model (w x , a x , and a c ), the initial 
value of the DE-DM ratio, Ri, should be such that 



Ri > i?_ 



(31) 



otherwise the evolution of the DE-DM ratio is not de- 
scribed by Eq. (|2T|) . A similar constraint exists for the 
late time value R^: 



Rno > Rn 



(32) 



to ensure that the value Ro is included in the range of 
values allowed by the exact solution (|27| . 

Further constraints on the values of the interaction pa- 
rameters can be obtained from Eqs. (|3"Tj) and (f3"2"| . Ac- 
cording to Eq. (JTTJ) , if the DE contribution is to be small 
at early times, we can write 



5 



and then, for a given value of Ri, Eq. (|3Tj) becomes the 
constraint 



Ri 



1 + R, 



■(1 



(34) 



By using Eq. (|18p . i.e. under the assumption that the 
DM contribution is small at late times, we find the com- 
panion expression of Eq. ()33|1 , which is 



R n 



Vb 



2Ci 



and then Eq. ([5"!?]) becomes the constraint 
a c < 1 - (1 + i?o) a x . 



(35) 



(36) 



If we want an upper limit on Roo, i.e. Roo < Roo.m, 
where Roo,m is some maximum value, then a new con- 
straint arises from Eq. (|33|) . 



(37) 



a c > 1 - (1 + R 

oo,m 



The inclusion of the constraint equations (|34|) . (|36j) . 
and (f3"T|) in Fig. Q] indicates that only a small region of the 
parameter space may be compatible with observations. 
The examples in Fig. Q] correspond to Ri = 0.1, i?o = 3.4 
and i?oo,m = 10. 

It can be verified that the above results are in agree- 
ment with the results presented in [J, [|| 0, H, l3ll l32j . 
In particular, the diverse cases presented in those papers 
are explained in a unified way by the exact solution (|27p . 
and our approach provides very simple expressions for the 
analysis of the parameter space. For example, it readily 
explains the troublesome features encountered in Figs. 7, 
8 and 9 of [|. 

In Fig. [2] we show examples of the evolution of the DE- 
DM ratio for fixed values Roo, Ro, and w x . The curves 
correspond to different values of a x , and the values of 
a c were determined from Eq. (|35p . Finally, we show in 



Fig.[3]a typical example of a phase space with parameters 
in the allowed region of Fig. [TJ 

We are assuming that the initial conditions are set at 
the onset of matter domination, and that the DM to 
baryonic matter ratio is the same as in the standard case, 



Pc 



(38) 



Note that the above is just an approximation, since the 
DM component in the interacting case does not evolve 
exactly as a -3 , see Eq. (|29|) . Then the initial conditions 
of the DM and DE contributions, with the help of the 
Friedmann constraint ([9]), are determined from 



6 

%i + -=Vi 
5 



1 



(39) 



Some examples of numerical solutions of the equations 
of motion are shown in Fig. [3] for different initial condi- 
tions. All of them represent valid solutions that start at 



z 
of 



0.1 



0.01 



0.001 



0.0001 




oL = 0.0850 

= 0.( 

= 0.( 
ACDM 
R = 3.4 




N 



10 



FIG. 2: Evolution of the DE-DM ratio R(N) according to 
Eqs. (|27[) and (|16[) . The chosen values of the various param- 
eters are Ro — 3.4, Rao = 10, whereas the value of a c was 
determined from Eq. (|35[) for the given values of a x . For com- 
parison, the green (dashed-dotted) line represents the stan- 
dard ACDM case. The dotted line corresponds to exactly 
etc = 0, for which case R-oo = and the earlier evolution is 
very similar to that of standard ACDM. 



a matter dominated epoch and end in a final state with 
finite DE-DM ratio. This final state is uniquely deter- 
mined by the values of the free parameters of the model. 

The phase space portrait reveals some other aspects 
of the interacting model under consideration, apart from 
the critical points and their stability studied above. It 
shows the line R-oo = x/y, which approximates very 
well the heteroclinic curve that connects the (unstable) 
critical point at the origin (actually, this is the baryonic 
dominated critical point discussed in Sec. lIII A|) with the 
(saddle) point A. It is then apparent that any curve with 
initial DE-DM ratio Ri < R-oo leads to trajectories in 
which the DE component is negative and the DM grows 
without bound. 

Therefore, the constraint equation (131|) is not only nec- 
essary for the trajectory to be described in terms of the 
DE-DM ratio in Eq. (|27|) . but it is also needed to have 
a reasonable evolution of the universe after the onset of 
matter domination. 

The initial value Ri is related to the e-fold number Ni 

by 

Ri = - — + — tanh [3C 2 (Ni - NJ] , (40) 
a x ol x 



where Aq is determined from Eq. (|28|) . so that there is 
a one-to-one correspondence between Ri and Aj. Thus 
different initial conditions result in different times for the 
appearance of a matter dominated epoch. Different ini- 
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FIG. 3: The phase space of Model I in the case 5t x = 0.15, and 
a c — 0.05. The (red) circles are the critical points A and B 
(see Table HI, and the green dashed line is the (heteroclinic) 
constraint x + y = 1 that connects them. The green dot- 
dashed line that connects the saddle point A to the unstable 
critical point at the origin is a good approximation to the 
heteroclinic trajectory between the two points. The red long- 
dashed line represents the ratio x/y — 0.1. The blue dotted 
line is the approximate Friedmann constraint (|39[1 . The black 
(solid) lines are numerical solutions of the equations of motion 
for different initial conditions. The trajectories with initial 
conditions on the right of the heteroclinic line end up at point 
B; if the initial conditions are on the left, the DE component 
becomes negative and the DM one grows without bound. 



tial conditions can have dramatic effects on the past his- 
tory of the universe, even if the final state can be ar- 
ranged to be the same in all cases. 



C. Special cases 

Our results apply for some special cases alreadypresent 
in the literature. For instance, if a x = <S C [5], then 
Eqs. (HHJ) suggest that a c < 1/4 for the critical points 
to be real (see also Eq. ([25)) ). Moreover, Eq. (fTT)) implies 
that & c -C 1 if proper matter domination is to appear 
and then R A < 1. Since Eqs. (H]) and together im- 
ply that Rb = 1/Ra, this model has difficulty to achieve 
an asymptotically constant DE-DM ratio relevant to the 
coincidence problem. 

Another example is a x = |9], for which the critical 
points, from Eq. (|15[) . are xa — a c and xb = 1- However, 
Fig. [1] indicates that any reasonable model necessarily 
needs a x ^ to alleviate the coincidence problem. 

For completeness, we can also have the case a c = 0. 
Then there is proper matter domination for any value 
< 6t x < 1, because x a = 0, and the case is free from 



the problems related to a finite value of R-oo , see Figs. [T] 
[21 and [3] The coincidence problem can be addressed 
if the only non-zero interaction parameter is given an 
appropriately small value, as Rb — 1/5^. 



IV. MODEL II. Q = 3 (T x p x + F c p c ) 

In the simplest model of the reheating process after in- 
flation in the early universe, one assumes that the oscil- 
lating inflaton field (j> behaves like a matter fluid that de- 
cays into relativistic particles; the decay is parametrized 
by a constant decay width [25]. In our notation, 
Q = T^ptf). Motivated by this, and by similar models 
for curvaton decay [2(| and for decay of DM to radia- 
tion we arrive at Model II: Q — 3(T x p x + T c p c ), 
where the F/ are constant decay widths. Unlike Model I, 
the Q here is not constructed a priori for mathemati- 
cal simplicity, and the dynamics are considerably more 
complicated as a result. 

The variable z becomes 



(41) 



which typically grows in an expanding universe, and di- 
verges in the limit H — > 0. Because of this, it is conve- 
nient to define the new variable 



H + H ' 



(42) 



which allows us to compactify the evolution of z. Here, 
Hq denotes the current value of the Hubble parameter. 
Early times correspond to u — > 0, and late times to u — > 
const (if H < 0, then u — ► 1). As in Model I, we redefine 
the interaction variables as 



ai{u) :-- 



7/:=Tf- (43) 



1 - u (-w x ) 

Notice that Uq = 1/2 and then aj = 7//(— w x ). Thus 
3u (~f x x + j c y) 



1 



3[a x (u)x + a c (u)y] , (44) 



which is a time-dependent version of the Model I expres- 
sion HH). The equations of motion become 

x' — 3x (1 — x) — 3 [a x (u) x + a c (u) y] , (45a) 
y' = -3xy + 3[a x (u)x + a c (u)y] , (45b) 



-- — (1 + w x x)u(l - u) . 

ZW r 



(45c) 



where a prime again denotes d/d(—w x N). Note that the 
DE equation of state appears explicitly in Eq. (I45c[) be- 
cause the value of w x is necessary to know the time evo- 
lution of the cosmological model. 
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A. Critical points 

The system (1451) is autonomous, and, to begin with, 
Eq. (|45c|) admits the critical values u* = and it, = 1. 

If = 0, then Eqs. (|45a[) and (|45b[) suggest the criti- 
cal values: (a) a;* = 0, with to be determined from the 
Friedmann constraint (|9]); and (b) x* = 1 and = 0. 
These are the expected critical points from the general 
discussion in SecHH see Eqs. (fT0|) . The stability of these 
points can be established by standard methods, from 
which we find that (a) is unstable, while (b) is of sad- 
dle type for — 1 < w x < and stable for w x < — 1. 

On the other hand, for u* = 1 the only possibility, 
again from the general discussion in Sec. [Til is x*+y* = 1, 
and then the critical value x* is determined from Eq. (fl5|) 
but with variable ai. Equation (|15j) has solutions 



a; ± (u) = C\(u) =F C 2 (u) + a c (w) 



(46) 



where C/(tt) are defined as in Eq. (|16p . with a/ — > a/(u). 

For u — > 1 (that is, a/ — > oo), the asymptotic approx- 
imate solutions are 



2 2w x (l - u) 

± 1 ( |7x - 7c 

2 



7a: ~7c 



[(Tz - 7c) ± |7a - 7c 



7x +7c 

Ix - 7c 



(47) 



The asymptotic values in the above equation depend on 
the values of the interaction parameters. In the case (7^ — 
7 C ) > 0, we find 



—00 
7c/ (7c -7a) 



and for (-/ x — j c ) < 0, we find 
7c/ (7c ~7x) 



vt 



.</* 



l-xl 



i-xi 



(48) 



(49) 



The critical points (xf.y+.l) are always unstable, 
whereas (x~,y~, 1) are always stable. The finite critical 
points of the system (|45() and their stability are summa- 
rized in Table HU 

It should be stressed that Model II is a time-dependent 
generalization of Model I, and then there should be good 
agreement in the formulas of the two cases. For in- 
stance, a quick comparison between Eqs. (jTTJ) and (JTHJ) 
and Eqs. (|48p and (|4"9"| . quickly shows this is the case. 

Point A represents matter (DM and baryons) domina- 
tion at early times, with no contribution from DE (unlike 
in Model I). Point C is a late-time attractor only in the 
case ("f x — 7 C ) > 0, and then we also require 7 C < if the 
energy density of the dark fluids is to be positive at the 
critical point. Furthermore, C is a scaling point, since 
the asymptotic DE-to-DM ratio is 



R 



x 
y 



7c 

Jx 



r x 



(50) 



which only depends on the ratio of the decay widths. This 
asymptotic ratio also shows that the interaction rates T/ 
should have opposite signs if both dark densities are to 
remain positive at late times. Thus 



> and r c < , 



(51) 



are necessary conditions to have a finite and positive late- 
time attractor in the model. 

Apart from finding the critical points, we need to know 
the behavior of the system at early times, since Model II 
can also exhibit a negative DE component, as shown in Q 
in the case T x = 0. 

In principle, we would need to scan exhaustively the 
3-dimensional phase space. A short-cut we will take is to 
find the points x at which x'(x) = 0, see Eq. (|45a|) . for 
fixed values of the variables y and u. The result is 



x±(u) 



1 r 



(1 - a x ) ± V(i -a x ) 2 -4a c y 



(52) 



In the early universe, u — > 0, one possibility is x + = 1, 
but the interesting solution is 



a c y- 



(53) 



Clearly, the value of x above marks the point at which 
the x-component of the phase space velocity x' changes 
sign. 

If the evolution of the cosmological system departs 
from an early unstable point corresponding to non- 
negative dark components, then we must impose the con- 
dition x~ > 0, otherwise the DE variable will approach 
point A from below, x — > 0" . This generalizes the r x = 
result of [7| to any value of T x and w x (see also the dis- 
cussion below). As a consequence, a non- negative DM 
interaction T c > is necessary to ensure a positive DE 
density at early times. 

It is now clear that we cannot, in general, find a version 
of Model II in which negative values of the dark energy 
densities are consistently avoided at both early and late 
times. In this sense, Model II can only be used if we 
restrict it to apply either (a) from the beginning only up 
to the present time, or (b) from some finite time (e.g. 
recombination) onwards. 



B. Special cases 

The first special case is = T c . Following the dis- 
cussion in the previous section, it is necessary to have 
positive interaction rates to have a proper early matter 
era. For the late Universe, the critical points are 



1 , 1 



= o ± o\/l-4a c (") 



(54) 



and then a negative value of T c is required for the critical 
points to be real. But even in this case the values of x^ 
are not finite in the limit u — * 1; hence, the simple case 
T x = T c cannot provide a realistic model. 
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Point 




y* 




Existence 


Eigenvalues 


Stability 


A 





y* 





all 7 C , J x 


0; 3; -3/(2w») 


Unstable if w x < 


B 


1 








all 7c , 7^ 


-3(l + tB a! )/(2io«); 
-3; -3 


Saddle if -1< w x < 
Stable if w x < — 1 


C 


-7c/ (7* ^Tc) 


7*/(7z - 7c) 


1 


7c + Ix 


3 7a! — 7c (l+Ulat) . 
2w x 7a, — 7 C 5 

37c/(7* -7c); 
—00 ■ sgn(7o, - 7c) 


Stable if — tiijic < 1, and xc > 0, and j x > 7c 
Unstable otherwise 



TABLE II: Critical points of Model II. The early universe is u —* 0, the late universe is u — » 1, and 71 = Vi/Hq. The 
co-eigenvalue for point C appears due to the limit u — > 1. As discussed in the text, see for instance Eq. (|5ip . 7^ > and 7 C < 
are required to have non-negative dark components at late times; those same conditions directly imply that point C is stable. 



In our numerical experiments we have found that the 
case ("f x — 7 C ) < has interesting properties. Firstly, a 
problematic early evolution can be avoided if we choose 
7c > 0; note that 7 X may positive or negative, as long as 
Ix < 7c- As u — ► 1, the point C is unstable and there are 
two possibilities for its late time evolution, see Eq. P5|) . 

The first one arises if at some time the DE variable x 
is to the right of point C, so that the system can freely 
follow the late-time attractor x~ — 00. This case is not 
interesting, as the Friedmann constraint demands then 
that the DM variable y — > — 00. 

The opposite case corresponds to the DE variable lo- 
cated to the left of point C, in which case x — > —00 and, 
because of the Friedmann constraint again, the DM vari- 
able grows without bound y — > 00. 

Whether the solution of Eqs. ([45]) is to the right or to 
the left of the unstable point C depends on the initial 
conditions and the values of the interaction parameters 
Ti. For instance, if the universe were described by this 
variant of Model II, then we would be to the right (left) 
of point C if xo > x+ (xq < x+ ). 

The special case previously considered in Q corre- 
sponds to = 0, and then it is a simple realization 
of the case just described above. According to our dis- 
cussion above, the model needs T c > in order to have 
a positive early DE component, in agreement with Q. 
(Note that our interaction parameters are defined with 
an opposite sign to those of [7|.) 

If r c < 0, the early evolution of the Universe is prob- 
lematic, but there exists the late time attractor x —* 1 
and y — > 0, which is precisely the case in Q. However, 
the case F c > is not a better option, because point C 
then becomes an unstable critical point for which x* = 1 
and y* = 0. According to its present conditions, our uni- 
verse would be at the left of point C and its late time 
evolution would result in x — > — 00 and y — > 00. 

Actually, Eq. (|lb[) has an exact solution if T x = 0, 

Pc = Pcoa~ 3 exp [3T c (t - f )] , (55) 



which tells us that the DM energy density is always posi- 
tivie, regardless the sign of r c . It also confirms the expec- 
tations discussed above: (a) if r c < 0, the DM compo- 
nent may scale at a rate faster than the usual a -3 at early 
times, but it exponentially vanishes at late times; (b) if 
r c > 0, the DM component may scale at a rate slower 
than the usual a~ 3 at early times, but it exponentially 
grows at late times. 

As for the special case T c = 0, we see that the dark 
components are well behaved, because both cases T x > 
and T x < lead the Universe to the (unstable) critical 
point 2+ = at early times. 

The differences appear at late times. If T x < 0, the 
discussion about the general case (7^ — 7 C ) < also ap- 
plies here. For example, our Universe would be presently 
located at the right of point C, and then its final fate 
would be x — > 00 and y — > — 00. 

If Tie > 0, then point C is stable but represents a DM 
dominated stage. If this were the case of the present Uni- 
verse, this would mean that the present DE dominated 
epoch is a transient phenomenon and that the Universe 
would eventually be dominated by DM again. 

The above statements can be clearly seen from the ex- 
act solution of Eq. (fTc)) in the case T c = 0, 

p x = p x0 a- 3(1+w ^ exp [-3T x (t - to)] , (56) 

which is the companion solution of Eq. (|55[) . The tran- 
sient character of the DE dominated epoch is controlled 
by the interaction T^; a long enough DE era requires an 
appropriately small DE interaction. 

C. DE-DM ratio 

The evolution of the DE-DM ratio is governed by 
Eq. (f25|) . with 07 — » aj(u). We no longer have an exact 
solution because of the explicit dependence on u; how- 
ever, there are some semi-analytical results that can help 
us to understand the dynamics of the model. 
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The troublesome features of negative DE at early times 
can also be derived via the equation for R: 



-1,Yx = 0-2 



r = -3(w x h + r c + r x )R - 3r c - z 1x r 2 



(57) 



In the matter dominated era, we have \R\ <C 1 and <C 
H, so that p c oc a -3 and then H = 2/3i. Equation (157)) 
becomes 



I? 



2u> T 



with solution 



— 



3r r 



1 + 2w, 



-R - 3r c 



-t + cr 



(58) 



(59) 



where C is an integration constant that has to be chosen 
to impose appropriately the condition R — > as t — > 0. 
If w x < —1/2, the r c mode dominates over the C mode 
as t — » 0, and then the DE becomes negative for r c < 0. 

There is another constraint we may impose on the free 
parameters of Model II. If we require the DE-DM ratio to 
be a growing function at the present time, i.e., R'\o > 0, 
then according to Eq. (|25[) . we need 



7c > -%Rq 



Rt\ 



1 + Ro 



(60) 



where 5/(1/2) = 7/, see Eq. (1T3"|) . 

We show in Figs. |4] and [5] some numerical examples 
for the evolution of the DE-DM ratio obtained from 
Eqs. (|45|) . The different cases confirm the results on 
Model II as discussed above. 

In particular, if we require positivity of both dark com- 
poinents as t — > 00 and as t — > 0, then we require T c = 
and a positive small value of r^; the smaller T x is, the 
larger the maximum value reached by the DE-DM ra- 
tio before the Universe enters again into a late-time DM 
dominated epoch. 



V. CONCLUSIONS 

We have made a careful analysis of two simple mod- 
els of interaction between DM and DE. The first model, 
Model I, has an interaction term proportional to the Hub- 
ble parameter times a linear combination of the dark 
sector densities, and is one of the most studied in the 
literature. 

We developed different mathematical techniques to in- 
vestigate the properties of the model under simple but 
general enough assumptions. Our results recover those 
of previous studies, and we found new analytic expres- 
sions that clarify the limitations of Model I. 

To begin with, we absorbed the (constant) DE equa- 
tion of state w x into the equations of motion, so that the 
parameter-space is truly 2-dimensional; this very much 
simplified the study of the allowed values of the inter- 
action parameters (compare our Fig. [T] with the corre- 
sponding 3-dimensional figures in [j]). 
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FIG. 4: Examples of the DE-DM ratio R as obtained from 
the numerical solutions of Eqs. (I45|l . under the condition that 
all cases have xo ~ 0.7 and yo = 0.24 at N = 0. We took a 
fixed value of 7^ = 0.2 for the case w x = —1, and the values of 
7 C are as indicated on the plot. Negative (positive) values of 
DE at early times appear for negative (positive) values of j c 
(see the inset), and a finite late time attractor appears only 
if the condition (7^ — 7 C ) > holds. The values for which 
(7a; — 7c) < lead to x — » -co, and may also break the 
constraint (l60l). 




FIG. 5: The same as in Fig. [4] but now for the fixed value 
7c = 0.2. As expected, early evolution is positive because 
7c > (see the inset), but point C is unstable in all cases. 
The DE-DM ratio first diverges for xq < xt , but asymptotes 
to R — > ± — 1 at late times whatever the case xo > xt or 
x < xt, see Sec. IIV Bl and Eqs. (|49)l . 



The absorption of the equation of state is not a mere 
mathematical trick. The original equations of motion 
would have been 



x' = —3w x x (1 — x) 
y' = 3w x xy + z . 



(61a) 
(61b) 



We can see that Eqs. ([S]) are recovered by setting w x = 
— 1 in the above equations. Thus, the cases discussed are 
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in some sense isomorphic to the cosmological constant 
interacting case. 

This fact also shows the degeneracy between interact- 
ing models with a constant DE equation of state; if there 
is a successful model with a cosmological constant, one 
can find another one with a different value of w x . How- 
ever, the models are distinguishable in their particular 
evolution, as the DE equation of state should appear ex- 
plicitly in the final expression of the DE-DM ratio, see 
Eq. ([77]) , and other quantities. 

Even though our general assumption was w x < 0, it 
is clear that in practice we have to restrict ourselves to 
values of the DE equation of state that allow an accel- 
erating expansion at the present time, i.e. wtot < —1/3. 
It is also possible to allow phantom values w x < — 1, al- 
though it is not clear whether any physically consistent 
models exist in this regime. A comparison with obser- 
vations would then be required to distinguish among the 
different possibilities, but this is beyond the purposes of 
the present study. 

We presented for the first time the surprisingly simple 
equation (TT2")) for the DE-DM ratio, and its simple ana- 
lytic solution (|27p . We stress that the exact solution [)27p 
allows us to easily and directly uncover the limitations 
of this interaction model. The exact solution has only 
one free integration constant, and thus we need a careful 
choice of initial conditions ensure an evolution that re- 
mains close to that of the standard non-interacting cos- 
mology. 

Surprisingly restrictive limitations are needed on the 
interaction parameters, a c and a x . Our simplified study 
of the parameter space showed that the parameters 
should take very small values. This is also in agree- 
ment with other studies in which one can find a careful 
comparison of the model with cosmological observations 
(see for instance 0, i, [H, [H H^). We showed that 
the model that seems to work better is the simple case 
a c = 0. This model is almost indistinguishable from the 
standard ACDM at early times, but provides a finite DE- 
DM ratio at late times so that the coincidence problem 
is alleviated. 

Other interesting properties arise from the interacting 
Model II of Sec. IIVI Even though it can be thought of 
as a time-dependent version of Model I, the properties of 
the critical points differ significantly. 

First of all, its vanishing interaction variables at early 



times allow Model II to have a true matter dominated 
epoch, which in fact corresponds to a (unstable) critical 
point of the dynamical system. However, a positive DM 
interaction T c > is necessary in order to keep the DE 
component positive at early times. 

Also, the DE-DM ratio can be finite and positive at 
late times, but this requires the dark interactions to have 
opposite signs and to comply with the condition T x > T c . 
It is then apparent that we cannot, in general, have a 
realization of Model II in which the dark components are 
well behaved for all times and the coincidence problem is 
addressed. 

There is though a simple model that may be realistic 
and corresponds to T c — 0. This model can satisfy all the 
constraints and the DE interaction can be adjusted so 
that the Universe can have appropriate matter and DE 
eras. The only difference is that DE domination is a 
transient event: the Universe would eventually go back 
to a DM dominated epoch at late times. 

Finally, we note that the problem of negative dark sec- 
tor densities is not the only problem with Models I and 
II: the curvature perturbation also has a non-adiabatic 
instability on large scales @ . As pointed out in [?|] , both 
of these pathologies in the interaction models are related 
to the way in which we treat DE, i.e. as a fluid with con- 
stant w x . If w x is allowed to vary in the early universe 
(as happens with scalar field DE), then the pathologies 
can be avoided. 



Acknowledgments 

GC-C is supported by the Programme Alban (the Eu- 
ropean Union Programme of High Level Scholarships for 
Latin America), scholarship No. E06D103604MX, and 
the Mexican National Council for Science and Technol- 
ogy, CONACYT, scholarship No. 192680. The work of 
RM was supported by the UK's Science & Technology Fa- 
cilities Council, and by a Royal Society exchange grant. 
The work of LAU-L is supported by CONACYT (56946), 
DINPO and PROMEP-UGTO-CA-3. RM thanks the 
Departmento de Fi'sica at the Universidad de Guana- 
juato for hospitality during part of this work. LAU-L 
would like to thank Diana Juarez and Daniel de la Torre 
for very helpful discussions. 



[1] J. Dunkley et al. (WMAP) (2008), 0803.0586. 

[2] M. Tegmark et al. (SDSS), Phys. Rev. D74, 123507 

(2006), astro-ph/0608632. 
[3] W. J. Percival et al., Mon. Not. Roy. Astron. Soc. 381, 

1053 (2007), 0705.3323. 
[4] M. Quartin, M. O. Calvao, S. E. Joras, R. R. R. Reis, 

and I. Waga, JCAP 0805, 007 (2008), 0802.0546. 
[5] L. P. Chimento, A. S. Jakubi, D. Pavon, and W. Zimdahl, 

Phys. Rev. D67, 083513 (2003), astro-ph/0303145. 



[6] A. P. Billyard and A. A. Coley, Phys. Rev. D61, 083503 

(2000), astro-ph/9908224. 
[7] J. Valiviita, E. Majerotto, and R. Maartens, JCAP 0807, 

020 (2008), 0804.0232. 
[8] C. G. Boehmer, G. Caldera-Cabral, R. Lazkoz, and 

R. Maartens, Phys. Rev. D78, 023505 (2008), 0801.1565. 
[9] Z.-K. Guo, N. Ohta, and S. Tsujikawa, Phys. Rev. D76, 

023508 (2007), astro-ph/0702015. 
[10] C. Wetterich, Astron. Astrophys. 301, 321 (1995), hep- 



11 



th/9408025. 

[11] L. Amendola, Phys. Rev. D60, 043501 (1999), astro- 
ph/9904120. 

[12] W. Zimdahl and D. Pavon, Phys. Lett. B521, 133 (2001), 

astro-ph/0105479. 
[13] G. R. Farrar and P. J. E. Peebles, Astrophys. J. 604, 1 

(2004), astro-ph/0307316. 
[14] R.-G. Cai and A. Wang, JCAP 0503, 002 (2005), hep- 

th/0411025. 

[15] G. Olivares, F. Atrio-Barandela, and D. Pavon, Phys. 

Rev. D71, 063523 (2005), astro-ph/0503242. 
[16] H. M. Sadjadi and M. Alimohammadi, Phys. Rev. D74, 

103007 (2006), gr-qc/0610080. 
[17] C. Quercellini, M. Bruni, A. Balbi, and D. Pietrobon 

(2008), 0803.1976. 
[18] J.-H. He, B. Wang, and E. Abdalla (2008), 0807.3471. 
[19] X.-m. Chen, Y. Gong, and E. N. Saridakis (2008), 

0812.1117. 

[20] N. Dalai, K. Abazajian, E. E. Jenkins, and A. V. 
Manohar, Phys. Rev. Lett. 87, 141302 (2001), astro- 
ph/0105317. 

[21] E. Majerotto, D. Sapone, and L. Amendola (2004), astro- 



ph/0410543. 

[22] H. Garcia-Compean, G. Garcia- Jimenez, O. Obregon, 
and C. Ramirez, JCAP 0807, 016 (2008), 0710.4283. 

[23] M. R. Setare and E. C. Vagenas (2007), 0704.2070. 

[24] N. J. Nunes and D. F. Mota, Mon. Not. Roy. Astron. 
Soc. 368, 751 (2006), astro-ph/0409481. 

[25] M. S. Turner, Phys. Rev. D28, 1243 (1983). 

[26] K. A. Malik, D. Wands, and C. Ungarelli, Phys. Rev. 
D67, 063516 (2003), astro-ph/0211602. 

[27] R. Cen (2000), astro- ph/0005206. 

[28] E. J. Copeland, A. R. Liddle, and D. Wands, Phys. Rev. 

D57, 4686 (1998), gr-qc/9711068. 
[29] S. del Campo, R. Herrera, and D. Pavon, Phys. Rev. 

D78, 021302 (2008), 0806.2116. 
[30] S. del Campo, R. Herrera, and D. Pavon, JCAP 0901, 

020 (2009), 0812.2210. 
[31] J. D. Barrow and T. Clifton, Phys. Rev. D73, 103520 

(2006), gr-qc/0604063. 
[32] G. Olivares, F. Atrio-Barandela, and D. Pavon, Phys. 

Rev. D77, 063513 (2008), 0706.3860. 



